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ABSTRACT 

We present a method from an X-ray observation of a galaxy cluster to measure 
the radial profile of the dark matter velocity dispersion, o"dm, and to compare 
the dark matter "temperature" defined as fJ^rn^a'^^/kB with the gas temperature. 
The method is applied to the XMM-Newton observation of Abell 1795. The 
ratio between the specific energy of the dark matter and that of the intracluster 
medium (ICM), which can be defined as /3dm in analogy with Pspec, is found to 
be less than unity everywhere ranging ~ 0.3 — 0.8. In other words, the ICM 
temperature is higher than the dark matter "temperature", even in the central 
region where the radiative cooling time is short. A /3dm value smaller than unity 
can most naturally be explained by heating of the ICM. The excess energy of 
ICM is estimated to be ~ 1 — 3 keV per particle. 

Subject headings: dark matter — X-rays: galaxies: clusters 
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1. Introduction 

Early X-ray imaging observations with the Einstein observatory and ROSAT showed 
that in the central regions of clusters of galaxies the radiative cooling time is shorter than 
the age of the universe (e.g Canizares, Stewart, & Fabian 1983). As a result, the intra- 
cluster medium (ICM) should cool down to form a cold {T < lO^K) gas phase inducing a 
global inflow of gas. This "coohng flow" (see Fabian 1994 for a review) picture has been 
extensively discussed and formed a basic assumption in many arguments. The low resolution 
spectroscopy in 0.5-2 keV by ROSAT showed that in some clusters the ICM temperature 
actually decreases towards the center (e.g. Bohringer et al. 1994; David et al. 1994; Allen 
& Fabian 1994). Higher resolution spectroscopy in 0.5-10 keV with ASCA, however, can 
not be fully understood with the conventional cooling flow model. ASCA spectra of cooling 
flow clusters can be well explained by a two (hot and cool) phase plasma without significant 
excess absorption features (e.g. Ikebe et al. 1999; Makishima et al. 2001). A naive cooling 
flow model predicting a range of temperatures with intrinsic absorption could also fit the 
ASCA data but generally produce worse chi-square results (e.g. Allen et al. 2001). Most 
recently, very high resolution spectroscopy with XMM-Newton/RGS unambiguously show 
that there is very little X-ray emission from gas cooler than certain lower cut-off tempera- 
tures of ~ 1 — 3 keV (Tamura et al. 2001; Kaastra et al. 2001; Peterson et al. 2001). Unless 
a large amount of cooled gas or the metals in the cold gas are hidden (Fabian et al. 2001a), 
there must exist a heating mechanism that prevents the ICM from radiative cooling. 

The necessity of global heat input into the ICM in addition to gravitational heating has 
also been pointed out from the break of the self-similarity between dark matter and ICM, 
which is most clearly demonstrated in the X-ray luminosity-temperature relation {L — T 
relation). A simple scaling analysis (Kaiser 1986) suggests a relation of L oc T^, while 
observation shows LcxT^ (e.g. Edge & Stewart 1991; David et al. 1993; White et al. 1997; 
Wu et al. 1999). The break of the self-similarity is also seen in the entropy vs temperature 
relation. Ponman, Cannon, & Navarro (1999) showed that cooler systems (T < 4 keV) have 
entropies higher than achievable through gravitational collapse alone. Simulation studies 
show that the observed relations can be reproduced, if there is enough non- gravitational 
heat input into the ICM by feedback from galaxies (e.g. Metzler & Evrard 1994; Bower 
et al. 2001) or preheating before the cluster formation (e.g. Navarro et al. 1995; Tozzi & 
Norman 2000). 

In order to shed some new light onto these "cooling flow phenomena" and "break of self 
similarity", we, in the present paper, perform a comparison of the temperature distribution 
of the ICM to the distribution of the velocity dispersion of the dark matter. A parameter, 
Pspec = '^^ail {^'s.T / l^'TTip), is oftcu uscd as a measure of the average kinetic energy per unit 
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mass in galaxies relative to that in the ICM. From observations of many clusters, the mean 
Pspec is ~ 1 with large scatter (e.g. Wu et al. 1999), indicating that the energy equipartition 
between galaxies and ICM is roughly achieved on average. In analogy with (3spec, we can 
introduce /^dm = (^DM/i^BT/fJ^mp) for comparison between the mean kinetic energy of the 
dark matter and that of the ICM, and define the dark matter "temperature" as Tdm = 
Hm^a'^y^/k-Q- Note that the above definition of "temperature", by means of proton mass 
instead of the actual mass of dark matter particles, is only for the sake of comparison with 
the gas temperature. Therefore we put "temperature" in quotes. We obtain, in this paper, 
the radial profile of the /3dm value observationally for the first time. 

In Sect. 2, we describe the method of measuring the dark matter velocity dispersion in a 
cluster of galaxies from an X-ray observation. We applied the method to the XMM-Newton 
data of a prototypical cooling fiow cluster, Abell 1795 (hereafter A1795), which is located at 
z=0.0616 (Struble & Rood 1987). The X-ray data analysis and results are presented in Sect. 
3. Discussion and summary arc found in Sect. 4 and Sect. 5, respectively. Throughout the 
paper, the Hubble constant is given as 70 /i7o km s~^ Mpc~^, and a fiat universe {Vt^^ — 0.3, 
^A,o = 0.7) is assumed. At the redshift of A1795, 1 arcsec corresponds to 1.19 kpc. 



2. Method of measuring dark matter velocity dispersion 

We use a simphfied model cluster being composed of a hot plasma ICM with a temper- 
ature of 10^ — 10^ K, and dark matter made of collisionless particles. Under the assumptions 
of spherical symmetry and hydrostatic equilibrium, the ICM distribution is described by 

GM _ kBTg /dlnrig dlnTA 
~ ~ fimp \d\nR ^ dlnRJ ' 

where M(< R) is the total gravitating mass within a sphere of radius R, /cb is the Boltzmann 
constant, G is the Gravitational constant, ng{R) and Tg{R) is the density and temperature 
of the ICM, respectively. Prom an X-ray observation, ng{R) and Tg(i?) are measured and 
M(< R) can be obtained via Eq. (1). When the dark matter particles tracing the same 
gravitational field are in steady state, they obey the Jeans equation 

GM _ 2 f dhipDM dhia^\ , , 

'''''' \dh^ ^ ~dh^ J ' 

where (Jdm(-R) is the one-dimensional radial velocity dispersion, and is the mass density, 
which is given as pdm = j:^^ ~ jim^rig. Therefore, once the total gravitating mass, M, 
and the gas density, rzg, profiles are obtained from an X-ray observation, Eq. (2) contains 
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only one unknown parameter, (7dm(-R)- Equation (2) can be rewritten in the form of an 
ordinary differential equation for cr^j^ as 



which can be solved numerically under a given boundary condition to derive the velocity 
dispersion profile of the dark matter. 

Although a steady state dark matter distribution is assumed above, a dark matter 
density profile should actually be growing as matter is falling onto the system from outside. 
A numerical simulation by Pukushige & Makino (2001) shows that the dark matter halo 
grows in a self-similar way, keeping the density profile in the central region unchanged. This 
justifies the assumption at least in the central region. 

A similar technique as described above was used in determining the ICM temperature 
profile, using X-ray imaging data without spectroscopic information taken e.g. with the 
Einstein observatory. A model profile for the total mass is assumed and the temperature 
profile is derived from the equation of hydrostatic equilibrium so that the observed brightness 
profile is explained (e.g. Fabian et al. 1981; Hughes 1989). Here we solve the Jeans equation 
instead in order to derive the dark matter velocity dispersion. 



An ideal opportunity to study the mass profile as well as velocity dispersion profile 
of dark matter is provided by the XMM-Newton observation of a prototypical cooling fiow 
cluster, A1795. Early X-ray imaging observations with the Einstein observatory and ROSAT 
showed that the X-ray emission from A 1795 has almost circular symmetry with a small 
elongation along the north-south direction, which indicates that the cluster is well relaxed 
dynamically (Jones & Forman 1984; Boute & Tsai 1996). The radial brightness profile shows 
a huge central excess above a prediction from an isothermal /9-model profile, and the central 
excess luminosity gives a mass deposition rate of ~ 250h^Q Mq yr~^ based on the standard 
coohng fiow interpretation (Edge et al. 1992; Briel & Henry 1996). Prom the 0.5-10 keV 
spectrum taken with ASCA, however, a significantly smaller mass deposition rate is obtained 
(~ 66/iTO Mo yr"^ Fabian et al. 1994; ~ TO/i^ yr"^ Xu et al. 1998; ~ 150/1^ Mq 
yr~^ Allen et al. 2001). Instead of the cooling fiow model, the ASCA spectrum can be best 
described with a two temperature model without excess absorption (Xu et al. 1998). 

The XMM-Newton observation of A1795 was carried out during the performance- verification 
phase on June 26 2000 with total observing time of ~ 50 ksec. The ICM temperature pro- 




GM _ CTpM dppM 

R^ Pdm dR 



(3) 



3. 



A1795 and the XMM-Newton observation 
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file obtained with EPIC MOS as well as EPIC PN were already reported by Tamura et al. 
(2001) and by Arnaud et al. (2001), showing that it is almost isothermal at ~6 keV in the 
2-10 arcmin radius region, while it decreases towards the center reaching the minimum tem- 
perature at ~ 3 kcV. The Chandra observation gives consistent results (Ettori et al. 2002). 
A high resolution spectrum from the central region obtained with XMM-Newton/RGS has 
been analyzed by Tamura et al. (2001), showing a lack of signature of cool gas components 
below ~ 3 keV and an upper limit of 77 h^Q Mq yr~^ for the mass deposition rate is obtained 
if an isobaric cooling flow model is applied. 

In this paper, we analyze the data taken with the XMM-Newton/EPlC PN, MOSl and 
M0S2 to measure the total mass profile of A1795 and to derive the dark matter velocity 
dispersion profile using the method described in Sect. 2. 

4. Analysis of the XMM-Newton data and results 

4.1. Data screening and background subtraction 

A significant fraction of any XMM-Newton observation is often contaminated by the 
huge background count rate by soft proton fiares. In order to derive time intervals with 
stable background, we made the 0.5-10 keV band fightcurves individually with the PN, 
MOSl, and M0S2 data, and eliminated time periods where the count rate deviates from 
the mean value during quiescent periods by > 2o". The total usable exposure time thus left 
is 23 ksec, 31 ksec, and 34 ksec for PN, MOS 1 and MOS 2, respectively. In the present 
analysis, the single and double pixel events of the PN data and the single, double, triple, 
and quadruple events of the MOS data are used. 

The background consisting of the cosmic X-ray background and the high energy particle 
events was estimated from the Lockman Hole data taken during the revolution 7^ 70. In this 
sky region, the Galactic column density is as low as the A1795 field. After eliminating the 
contamination sources from the Lockman Hole data as well as the A1795 data, the Lockman 
Hole data are scaled so that the count rates of the background and the cluster data becomes 
the same in the 7-12 keV energy band, and in the radius range of 575-925 arcsec for PN and 
of 500-850 arcsec for MOS, where no significant X-ray signal from the cluster was detected. 
The background data thus obtained were subtracted from the cluster data. For PN, the 
out-of-time events (Striider et al. 2001) were also subtracted prior to the subtraction of the 
background. 
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4.2. Radial count rate profile 

Defining tiic center as tlie X-ray peak position, we have produced a projected radial 
brightness profile in the 0.8-10 keV band. The X-ray count rate was converted to the X- 
ray flux by correcting each X-ray photon with the instruments' responses that include the 
effective area and vignetting of the X-ray telescopes, and the quantum efficiency of the 
CCD cameras. Combining the data taken with PN, MOSl and M0S2, we thus derived 
the brightness profile, which is illustrated in Fig. 1. For comparison with the other X-ray 
observations, we quantified the radial profile using the /3-model profile given by 

-3/3+0.5 

(4) 

where Rc is a core radius and P is a. beta parameter. The (3 model has been often used to 
fit X-ray brightness profiles and generally gives a good representation for non cooling fiow 
clusters. Since the X-ray brightness profile of A1795 has long been known to show a central 
excess above a /9-model profile, we employed here a double /9-model profile that is the sum of 
two /3-model profiles to fit the 0.8-10 keV brightness profile within 848 arcsec. In the actual 
fitting, the model brightness profile is convolved with the point spread function of the X-ray 
telescope ^, whose half-power-radius is ~ 9 arcsec at 1.5 keV on-axis and significantly affects 
several of the central bins. Evaluating the goodness of the fit with a chi-square statistic, we 
obtained a good fit as shown in Fig. 1 with the best-fit parameters as summarized in table 
1. 



So 



R 
Rc 



4.3. ICM temperature and metallicity profile 

The temperature and metallicity profiles given by Tamura et al. (2001) were derived 
from the conventional annular spectral analysis, accumulating spectra from concentric an- 
nular regions in the projected 2-dimensional space. We derive here the temperature and 
metallicity profiles instead in 3-dimensional form from deprojected spectra. With each PN, 



where and a depend on 



^The PSF is approximated by an analytic function given as 

off-axis angle as well as X-ray energy (see Ghizzardi 2001). We assumed that the PSF is constant in the 

entire FOV and independent of the X-ray energy. As the representative, we employed the on-axis PSF at 
1.5 keV. Using the rc and a parameters of (rc[arcsec], a) = (5.37, 1.5), (4.72, 1.457), and (4.4545, 1.4035), 
for PN, MOSl, and M0S2, respectively, we obtained the average PSF, which is then used in the brightness 
profile fittings. 
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MOSl, and M0S2 data set, we first accumulate 13 annular spectra from the vignetting cor- 
rected count rate profiles in individual energy channels. The outer radii are 16, 32, 48, 64, 
96, 128, 192, 256, 384, 512, 640, 768, and 848 arcsec. Assuming spherical symmetry and 
that there is no cluster emission beyond 848 arcsec from the cluster center, we deprojected 
the annular spectra to spectra for 13 spherical shell regions whose radii correspond to the 
annular radii. Since no intrinsic absorption is found even from the central region (Tamura 
et al. 2001), the X-ray emission from the cluster is entirely optically thin, and we performed 
the deprojection process as follows. The relation between the energy spectra observed in the 
N annular regions, Sj{E) [j — 1 — N), and those emitted from the N shell regions, Ci{E) 
{i — 1 — N),\s given with a matrix relation as 

N 

= J]Mj,ia(£;) , (5) 

i=l 

where N — Vi m. the present case. The matrix elements are given as 

f [{rl, - rf )^/^ - (rf - rjfl' - (r^, - rf,,)^/^ + (rf - r-^,m/{rl, - rf) Hi > j, 

I if z < J, 

(6) 

where ri (rj) and ri+i (rj+i) are the inner and the outer radius of the i'th shell (j'th annular), 
respectively. Each element represents the volume fraction of the i'th shell contributing to 
the emission detected in the j'th annular. Since Mj^i is a triangular matrix, it can be easily 
inverted to M-~^^, and the deprojected shell region spectra can be derived from the annular 
spectra by 

TV 

Statistical errors in the individual energy bins are properly propagated through this opera- 
tion. 

Out of the 13 deprojected spectra thus derived from each PN, MOSl, and M0S2 data 
set, we used only the inner 10 deprojected spectra within a 512 arcsec radius in the analysis 
below so that the artificial cut-off of the X-ray emission outside the maximum radius of 
848 arcsec does not affect the resulting deprojected spectra within 512 arcsec. In addition, 
the signal to noise ratios of the spectra within 512 arcsec are high enough, and some 10% 
systematic error of the predicted background count rates introduce no significant systematic 
errors on temperature measurements. 

In order to derive the temperature and metallicity profile, we simultaneously fitted the 
PN, MOSl, and M0S2 spectra with a single-temperature plasma model in each deprojected 
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shell region. We employed the MEKAL model (Mewe et al. 1985, 1986; Kaastra 1992; 
Liedahl et al. 1995) and assumed the abundance ratio among different elements to have 
solar values. The fits are acceptable for all the ten deprojected shell regions. The depro- 
jected temperature and metallicity profiles thus derived arc illustrated in Fig. 2 and Fig. 3, 
respectively. Compared with the result from the conventional annular spectral analysis, the 
effect of the projection is clearly seen inside the ~60 arcsec radius. 

The temperature is consistent with being isothermal in the outskirts beyond ~200 kpc, 
while towards the center it decreases monotonically and reaches 2.8 keV in the innermost 
shell. Note that the central minimum temperature agrees with the cut-off temperature 
in a cooling flow model that describes the RGS spectrum (Tamura et al. 2001) and the 
temperature of the X-ray filament detected with Chandra (Fabian et al. 2001). For usage in 
the following sections, we obtained an analytical formula that approximates the temperature 
profile. We employed a formula given as 



T{R) ^To + Ti 



1 + 



R 



R, 



c,T 



-r; 



(8) 



which has been proposed by Allen, Schmidt, and Fabian (2001) to give a good description 
to radial temperature profiles in relaxed clusters. Leaving all the four parameters of Tq, Ti, 
Rc,T, and T) as free parameters, we fitted the temperature profile derived from the deprojected 
spectra. The parameters obtained from a minimization are To = 2.78lQ 3g (keV), Ti = 



Smt'oll (keV), i?c,T 



50.3^^r 



arcsec 



and 7] 



9 7Q+0.28 
^- '°-0.38 

1.98^0Q5, where the errors give the 90% 



confidence range for four parameters of interest (A^^ = 7.78). The best-fit function and the 
error ranges for the individual radii arc overlaid in Fig. 2, showing that the function gives a 
good representation of the temperature profile. 

In the metallicity profile (Fig. 3), the strong central concentration of metals found in the 
previous works (Tamura et al. 2001; Ettori et al. 2002) is confirmed. As for the temperature 
profile, the metallicity profile was also fitted with a similar function given as 



A{R) =Ao + Ai 



1 + 



R 



R, 



c,A 



-3/2 



(9) 



0.20;^-|!] (solar), Ai = 0.321°;}^ (solar), and i?c,A = 149l^^° (arcsec), and 



We derived Aq — w.^w_q 
the best-fit function is overlaid in Fig. 3. 
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4.4. Mass profile 

An immediate way to derive the mass profile is using Eq. (1) with the temperature 
profile measured in Sect. 4.3 and the density profiles of the ICM which can be obtained 
from the double /3-model profile fitted to the observed X-ray brightness profile in Sect. 
4.2. Assuming spherical symmetry, the projected count rate profile in the 0.8-10 keV band, 
E(0.8 — lOkeV), gives the deprojected radial emissivity profile, e(0.8 — lOkeV), which can be 
converted to the ICM density via the relation 

e(0.8 - lOkeV) = n^A(T, A; 0.8 - lOkeV) , (10) 

where A(T, A;0.8 — lOkeV) is the emissivity coefficient defined by temperature (T), metal- 
hcity {A), energy range (0.8-10 keV), and the instrument responses. 

Taking into account the error ranges in the temperature profile, we thus derived the 
total gravitating mass profile, which is shown in Fig. 4. The result roughly agrees with the 
mass profile given by Xu et al. (1998) derived from the ASCA data with the same method, 
although a shoulder-like structure found in their profile at ~ 100 kpc is less prominent in our 
profile. The difference in the mass profiles is mainly attributed to the different parameters 
of the double /?-model profile fitted to the brightness profile. 

4.5. Theoretical modeling to the mass profile 

The total mass profile thus derived can not be conveniently used in Eq. (3) for deter- 
mining the dark matter velocity dispersion, however. The term of dp^u/dR in Eq. (3) is 
calculated from the third derivative of the ICM density profile, i.e. the third derivative of 
the radial brightness profile which is approximated with the double-/? profile, and it gives 
an uneven profile at the radius where the two f3 profiles cross over. This causes an unstable 
behavior in the solution of Eq. (3) and introduces non-negligible systematic errors in the 
results. Moreover, the double-/? modeling on the brightness profile can only reproduce a 
mass profile with the flat core as long as the temperature profile is virtually flat near the 
center, and provides rather restrictive mass profile. 

Here we make an alternative approach from the derivation above. We start modeling 
the total mass profile with an analytical formula with a few free parameters that assures a 
smooth profile with Eq. (3). Combined with the temperature profile derived in Sect. 4.3, 
the mass profile is used to predict the X-ray surface brightness profile of the cluster to be 
compared with the data to adjust the free parameters. In the actual calculation, a given total 
mass profile, M(< R), and a temperature profile, T{R), is converted to the ICM density 
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profile from Eq. (1) as 



T(0) 



R 



GixmpM{< R) 
kT{R)R^ 



dR 



(11) 



where no is the central ICM density. The emissivity profile in a given energy range, e(-Ei, -E2), 
is then obtained by Eq. (10), where Wg, T, A are substituted by Eqs. (11), (8) and (9), 
respectively, which is then converted to the brightness profile, T^^Ei^E^). Being convolved 
with the PSF, the model brightness profile is fitted to the 0.8-10 keV count rate profile 
observed to determine the best-fit parameters in the mass profile model. Although the 
observational constraint on the ICM temperature profile is limited within 512 arcsec from 
the center, we apply Eq. (8) beyond the radius, where it describes a practically isothermal 
model. 

As for the mass profile, we apply two representative models having a fiat core and a 
cuspy core. The King model (King 1966) is a classic model for the gravitational potential 
structure characterized by a fiat core. We used an approximated formula of the King model, 
in which the density profile is given as 



1 



R_ 
Rc 



-3/2 



(12) 



while the integrated mass profile is given as 

R IR^ 



M{<R)^ ^TiRlPo 



In 



Rr 




R fR 



Rr V 



+ 1 



-1/2- 



(13) 



where po is the central density and Rc is the core radius (see e.g. Binney and Tremaine 1987). 
The predicted X-ray brightness profile from the King model is fitted to the 0.8-10 kcV count 
rate profile. Fitting parameters are the core radius (-Rc) and the central density (po) of the 
approximated King model, and the central ICM density (uq). The best-fit model and the fit 
residuals are shown in Fig. 5 and the parameters derived are summarized in table 2. The 
fit is not acceptable. In Fig. 4, the best-fit King model is compared with the mass profile 
obtained in Sect. 4.4, showing clearly that a central mass excess is necessary in additional 
to the King model mass. 

We then added another King-model component to account for the central mass excess 
to construct a double approximated King model, which still has a fiat core profile. Fitting 
parameters are now two sets with a core radius and a central density for each (it!c,i, Po,i, 
i?c,2, and Po,2), and the central ICM density (no). As shown in Fig. 5 and summarized in 
Table 2, a good fit was obtained. 
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Another analytical formula for the total mass profile we employed is the universal halo 
profile by Navarro et al. (1995, 1996, 1997; NFW model), which is characterized by a sharp 
central cusp. The density profile of the NFW model is given as 



where Rs is called scale radius. The integrated mass profile of the NFW model is given as 
(Makino et al. 1998; Suto et al. 1998) 



The X-ray count rate profile predicted from the NFW model is fitted to the data, where the 
free parameters are po, Rg, and uq. As illustrated in Fig. 5 and summarized in Table 2, the 
NFW model reproduces the data very well. 

For each best-fit mass profiles, the double approximated King model and the NFW 
model, the corresponding virial mass, Myir, and the virial radius, Rvh, are calculated and 
given in Table 2. The concentration parameters of the NFW model defined as c = Rg / Rvir is 
obtained to be 6.1. It is worth mentioning that the obtained c and M^[^ are consistent with the 
theoretically predicted c — Myir relation (e.g. Bullock et al. 2001). The double approximated 
King model and the NFW model obtained above are compared in Fig. 6, together with 
the mass profile obtained in Sect. 4.4. The three mass profiles are all consistent within the 
20-600 arcsec radius region, and they approximately follow M oc R^ '^ in this zone. In Fig. 
7, the double King model and the NFW model are compared in the mass density profile, 
together with the corresponding ICM density profile. As clearly seen in Fig. 6 and 7, the 
two theoretical model profiles differ most significantly in the very central region within ~20 
arcsec radius, where XMM can not well resolve the spatial structure. Ettori et al. (2002), 
using Chandra data, obtained the total mass profile in A 1795 with finer spatial resolution 
and found that the density approximately follows p oc in the central region. This 

profile is steeper than the King model (p oc it!°) and flatter than the NFW model (p oc R~^) 
near the center. Therefore, the two mass profiles give a conservative mass range. In the 
next subsection, for determining the dark matter velocity dispersion profile, we use both the 
double approximated King model and the NFW model for the total mass profile model, which 
would also be expected to give a conservative range of the velocity dispersion representing 
systematic uncertainty involved in the usage of a specific mass profile model. 




-2 



(14) 




(15) 
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4.6. Dcirk matter velocity dispersion profile 

We now calculate the dark matter velocity dispersion profile by solving the differential 
equation given in Eq. (3). As the total mass profile, M(< R), the double approximated 
King model or the NFW model profile obtained in Sect. 4.5 is employed, while the dark 
matter density profile is obtained from the total mass and the corresponding ICM density 
profile by Pdm(-R) — d^^^ — iim-pn^{R). Both M(< R) and Pdm(-R) are substituted in 

Eq. (3) and the differential equation is solved with the fourth-order Runge-Kutta method 
for various boundary conditions at the center, i.e. cr'^yii^)- 

Figure 8 shows a set of solutions of Eq. (3), when the best-fit NFW model is applied to 
M(< R). Note that there is only a weak observational constraint on the mass profile beyond 
~ 1 Mpc and the solution at the large radii depends on the validity of the assumed mass 
profile. Depending on the inner boundary condition, the dark matter velocity dispersion 
may fall to zero at small radii, or it may diverge to infinity. In order to select physically 
plausible solutions among them, we placed a conservative boundary condition that the ve- 
locity dispersion at the virial radius is greater than and less than the free fall velocity at 
this radius. 

The constraints on the velocity dispersion profile thus obtained are shown in Fig. 9 
individually for different mass models applied. The errors on the ICM temperature profile 
are properly taken into account to obtain the error bands. In the central region, there 
is a clear difference between the results from the double approximated King model and 
the NFW model as expected from their different behavior. The velocity dispersion is then 
converted to the "temperature" of the dark matter (/cTdm = (^dmA*"^?); which is compared 
with the ICM temperature in Fig. 10. The ICM temperature is greater than the dark matter 
"temperature" everywhere. Even in the central region where radiative cooling is expected to 
be most effective, the ICM temperature is significantly higher than that of the dark matter. 
The comparison of the temperatures can be more directly described by means of the /?dm 
value defined as jS^M = (y'^y^^m^ / kT = T^yi/T^ which ranges ~0.3-0.8 (Fig. 11). There is 
no sign that the ICM is cooled significantly below the dark matter "temperature" . In other 
words, the dark matter looks to form the temperature floor that limits the ICM temperature. 

As a matter of fact, what we derived by solving the Jeans equation is the velocity 
dispersion of the collisionless particles that includes the dark matter as well as the stellar 
component. As shown in Fig. 4, the stellar component makes a minor contribution to the 
total mass and the derived velocity dispersion profile shown in Fig. 10 is virtually that of the 
dark matter, except in the very central region, where the two components are comparable. 
If the galaxies alone are in a steady state, they should also obey the Jeans equation under 
the same gravitational potential. In Fig. 10, we overlay the velocity dispersion profile of 
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the member galaxies measured by den Hartog and Katgert (1996), and the stellar velocity 
dispersion of the cD galaxy provided by Blakeslee and Tonry (1992), after converted to 
"temperature" via /cb?"* = jim^al. They show roughly a consistent profile with that of the 
dark matter at least in the case when the total mass profile is given by the NFW model. 



5. Discussion 



5.1. Anisotropy of (Jdm 



Equation (2) that we have used to calculate the dark matter velocity dispersion profile is 
based on the assumption of isotropic motion of dark matter particles. However, simulation 
studies (e.g. Eke, Navarro, and Frenk 1998; Colin, Klypin, & Kravtsov 2000) indicate 
that the radial velocity dispersion, cr^, should be rather larger than the tangential velocity 
dispersion, erf = \{crl + (t|). The degree of anisotropy is often measured with 



a2 



A ^ 1 - -I , (16) 



and the Jeans equation is modified as 



GM ^_^,(d\^^dA^^ ^^^^ 



R ''''' \ dlnR dlnR 



Employing A = derived by Colin, Klypin, & Kravtsov (2000), which is 

as large as 0.5 at the virial radius, i?vir, and converging to at the center, we solved Eq. 17 
as done in Sect. 4.6. The o"dm profile thus derived for the case of the best fit NFW mass 
profile is shown in Fig. 12, overlaid with the A = solution given in Fig. 9. The cr^ value 
with the anisotropy is greater than that without the anisotropy by only ~6% at lOOkpc and 
at most ~40% around IMpc. 



5.2. Heating source 



The /3dm profile determined here from observations should provide information on the 
thermal history of the ICM. From numerical simulation studies, /9dm ~ 1 — 1.4 is expected, 
if there is no cooling or additional heating (e.g. Metzler and Evrard 1994; Navarro et al. 
1995; Bryan & Norman 1998). An obvious way to explain the /?dm value smaller than unity 
in A1795 is heating of the ICM. 
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As suggested from the break of the self-similarity between dark matter and ICM, there 
should have been non-gravitational heating acting globally. We, from our results, estimated 
the excess energy of the ICM over that of the dark matter particles as 

AEi< R) = ^{<kT> -iim^ < al^ >) , (18) 

where <> denotes mass weighted mean within radius, R. As shown in fig. 13, the excess 
energy, AE, thus derived is found to be ~ 1 — 3 keV per particle, which may be compared 
with theoretical model predictions. The amount of energy injection to the gas phase that 
explains e.g. the observed L — T relation depends on the period when the heating occurred. 
Heating prior to cluster collapse, "preheating" , needs 0.1 — 0.3 keV per particle (e.g. Navarro 
et al. 1995; Tozzi & Norman 2001), while heating after a cluster formation requires higher 
values of 1 — 3 keV per particle (e.g. Metzler & Evrard 1994; Loewenstein 2000; Wu, Fabian, 
& Nulsen 2000; Bower et al. 2001). Our results given above may indicate that the global 
non-gravitational heating that may cause the break of self similarity has happened mainly 
within a collapsed cluster. However, even if such non-gravitational heating that explains 
the global X-ray feature of the cluster is provided, the central region of A1795 has a short 
radiative cooling time (Fig. 14), and there must be another significant energy input to the 
ICM in the central region at the present epoch to prevent the ICM from cooling. 

As a possible energy source in the central region, we first consider gravitational energy 
of the member galaxies and stellar components therein. The kinetic energy of the random 
motion of stars can be partially transferred to the ICM by stellar mass loss. Gas supplied 
by stellar mass loss has velocities of the bulk motions of stars relative to the ICM, which is 
the sum of a galaxy motion and the motions of the stars in the galaxy, and is likely to be 
thermalized by interactions with the ambient gas. If the stellar component moving in the 
same gravitational potential has a similar velocity dispersion profile as that of the dark matter 
(Fig. 10), the gas temperature achieved from this process is expected to be comparable to the 
velocity dispersion of the dark matter. This process nicely accounts for X-ray emission from 
isolated X-ray compact elliptical galaxies (Matsushita 2001). The input rate of the kinetic 
energy of the gas from stellar mass loss may be simply estimated a.s E = l/2Msta,r^o'^ = 10^^ 
ergs/s, where M^tax is the total stellar mass of 1 x 10^^ M^, m=3 x 10^^ M© yr"^ (IO^^Mq)^^ 
is the stellar mass loss rate in unit time, and (t(=580±80 km/s) is the velocity dispersion of 
the galaxies from den Hartog & Katgert (1996; Fig. 10). (The stellar velocity dispersions in 
the galaxies are smaller and neglected here.) This is much smaller than the output energy 
in the central region by X-ray radiation in galaxy clusters like A1795, however. Thus for the 
case of galaxy clusters, we have to look for additional heat sources. 

Kinetic energy of the stellar component might be more efficiently transferred to ICM 
via magnetic fields. As pointed out by e.g. Makishima et al. (2001), the motion of stars 
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may amplify interstellar magnetic fields and reconnections of the fields may heat up the ICM 
rather efficiently. The galaxies must have lost their kinetic energies through interactions with 
the ICM and have gradually fallen inwards accumulating onto the central galaxy to form 
the cD galaxy. A deep optical image of the cD galaxy in A 1795 derived by Johnstone et al. 
(1991) shows a concentration of elliptical galaxies of various sizes and stars forming a largely 
extended envelope with 131 h^Q kpc effective radius, which strongly suggests the on-going 
formation process of the cD galaxy. Quantitatively, the total amount of dynamical energy 
of the stellar component in the member galaxies that has been lost in the past is estimated. 
The stellar component in the galaxies is assumed to have a negligible small potential (U) 
and kinetic energy (K) before the formation of the cluster, and the current energy of the 
stars is estimated to be fZ + i^ ~ —10®^ ergs. If the energy has been released over the last 10 
Gyr and has been spent in ICM heating, the heating luminosity is expected to be ~ 3 x 10'^^ 
ergs s^^. This amounts to the bolometric luminosity of the ICM within GOhj^ kpc (Fig. 14), 
and may be sufficient to sustain the thermal energy of the ICM against radiative cooling. 
This model predicts that the stellar velocity dispersion became smaller than that of ICM, 
i.e. Pspec < 1) which is consistent with the actual observed value in the central region (Fig. 
10). 

Alternatively, there may be sufficient non-gravitational heat input from an active galac- 
tic nucleus (AGN) powered by an accretion of low entropy gas at a cluster center (Ghurazov 
et al. 2002; Bohringer et al. 2002). Numerical simulations show that an outflow from an 
AGN may form hot bubbles of relativistic plasma rising with buoyancy, and the bubbles may 
uphft cold gas mixing with ambient ICM (Ghurazov et al. 2001; Quihs, Bower, & Balogh 
2001; Briiggen & Kaiser 2002; Basson & Alexander 2003). High resolution imaging obser- 
vations revealed ripples and shock features in the central region of Perseus cluster (Fabian 
ct al. 2003) and M87 (Forman et al. 2004), suggesting that the bubble energy may also 
be transfcrcd by sound wave to larger distances, and that shocks may be the major con- 
tribution of the energy dissipation. The heating mechanism is self regulated: the lower the 
entropy, the higher the accretion rate onto the central engine. A portion of the accretion 
power is dissipated back into the ICM to make its entropy high and regulate the accretion 
rate to achieve an equilibrium state. This process automatically prevents the persistence of 
cold and hence dense clouds. The cD galaxy of A1795 has a radio source, 4C26.42, and the 
existence of an AGN is evident. Using the physical state of the ICM in the center we can 
actually estimate the energy provided by the AGN by applying the Bondi accretion model. 
According to the well-known correlation of the black hole mass with the mass of the bulge 
component (Magorrian et al. 1998), the black hole mass is expected to be ~ 6 x 10^ Mq. 
Assuming that the gas profile is fiat in the center, we can use the measured values of ng=0.1 
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cm ^ and To—2.8 keV to obtain the Bondi mass accretion rate 

M = A7i0.25pooC~l{GMBHf (19) 

where poo and Cg^oo are the density and sound velocity outside the Bondi accretion radius. 
We find M ~ 0.4 x 10^'^ M,.; / yr . Under the standard assumption of 10% of the accretion 
energy to be dissipated, the output energy is found to be = O.lMc^ ~ 3 x 10^^ ergs/s. This 
amounts to the X-ray luminosity within the inner ~ 20 kpc region only. The Bondi accretion 
radius is estimated to be Rb ~ GM^u/Cs ^ 30 pc, much smaller than the resolution of the 
temperature and density structure that can be measured with XMM-Newton. If the ICM 
density is not uniform but is clumpy, the Bondi accretion rate should be significantly larger 
and the heating rate could also be larger than the above estimation. We can note that the 
ICM conditions might temporally vary and that we currently see a relatively low state. 

Note that there are other heating mechanisms discussed, which include e.g. a classical 
idea of thermal conduction from the hot outer regions (e.g. Takahara & Takahara 1979; 
Tucker & Rosner 1983; Bregman, & David 1988; Gaetz 1989), and the acoustic wave heating 
recently proposed by Fujita, Suzuki, & Wada (2004). Comparison of those model predictions 
with our results presented in the current paper would be very interesting. 

6. Summary 

We derived the dark matter velocity dispersion profile from an X-ray observation for 
the first time. Using the XMM-Newton EPIC data of A1795, we derived /3dm of ~ 0.3 - 
0.8, indicating that the ICM temperature is larger than the dark matter "temperature" 
everywhere. We also derived the excess energy in the ICM, AE, which is found to be 
~ 1 — 3 keV per particle. These observational quantity should be useful to provide a new 
clue to solve the heating mechanism in the cooling core and to understand the thermal 
history of a galaxy cluster. 

We thank the XMM team for providing the data and the data analysis tools. In par- 
ticular, we express our gratitude to Michael Freyberg for helping the data analysis. We also 
acknowledge FTOOLS. We thank Kuniaki Masai for helpful discussions. We are also grateful 
to Paul Lynum for critically reading this manuscript. 
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Fig. 1. — In the upper panel, the 0.8-10 keV X-ray count rate profile shown in crosses is 
fitted with a double /3-model (solid line). The dotted lines show the individual /3- model 
components. The fit residuals are shown in the lower panel. 
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Fig. 2. — The ICM temperature profile of A1795 derived from the spectral analysis 
with the XMM-Newton EPIC-MOS and EPIC-PN data. Results from the deprojected 
spectra are shown with diamonds, while the conventional annular spectral analysis gives 
the profile with crosses. The solid line shows the best-fit function given as T{R) = 



2.78 + 3.06 



1 + 



R{arcsec) 
503 



-1.98 



region with the function. 



(keV), while the gray hatched region shows 90% error 
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Fig. 3. — The metallicity profile of A1795. Results from the deprojected spectra are shown 
with diamonds, while the conventional annular spectral analysis gives the profile with crosses. 
The solid fine shows the best-fit function given as A{R) = 0.20 + 0.32(1 + {R/U9y)-^-^. 
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Fig. 4. — The grey hatched region illustrates the integrated radial profile of the total gravi- 
tating mass obtained via Eq. (1). The errors in the ICM temperature profile illustrated in 
Fig. 2 is properly propagated to obtain the error range at the individual radii. The vertical 
dotted line indicates the radius within which the ICM temperature measurements are avail- 
able, and isothermality is assumed beyond the radius. The dashed line shows the ICM mass 
profile calculated from integrating the ICM density profile, where the errors are smaller than 
the line width. The dot-dashed line shows the stellar mass profile, estimated from an J-band 
image obtained by Johnstone et al. (1991). We approximated the J-band surface brightness 

profile with a de Vaucouleurs law given as /i[mag arcsec"^] = 25.0 + 8.33 

and converted it to the mass profile assuming 5-1=2.28 and a mass-to-light ratio of 6.5 in 
i?-band. The solid line indicates an approximated King profile (Eq. 12 and 13) derived by 
fitting it to the data via Eqs. (10) and (11). 
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Fig. 5. — The upper panels show the 0.8-10 keV X-ray count rate profile (crosses) and the 
best-fit model profiles predicted from given total mass profiles. The mass profile assumed 
in the model for the left, middle and right panel is single approximated King model, double 
approximated King model and NFW model, respectively. The lower panels show residuals 
for the best-fit model profiles. 
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Fig. 6. — The scaled total mass profiles obtained in Sect. 4.4 (light gray band) and the two 
theoretical models, the double approximated King model (hatched region with oblique lines) 
and the NFW model (dark gray band) obtained in Sect. 4.5 are compared. The vertical axis 
shows the integrated mass multiplied by R^^ '^ . The vertical dotted fine is the same as in 
Fig. 4. 
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Fig. 7. — The upper curves show the mass density profiles of the total gravitating matter, 
while the lower curves show the ICM density profiles. The results with the double approxi- 
mated King model and those with the NFW model are illustrated with light gray band and 
dark gray band, respectively. The vertical dotted line is the same as in Fig. 4. 
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Fig. 8. — Solutions of Eq. (3) with different (Jdm(O) values, when the best-fit NFW model 
is used for the mass profile. Physically plausible solutions are illustrated in the solid lines, 
while the dashed lines are other possible solutions. The vertical dotted line is the same as 
in Fig. 4. 
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Fig. 9. — The velocity dispersion profiles of the dark matter are shown with the error bands. 
The gray hatched region and the hatched region with oblique lines are derived from the double 
approximated King model and the NFW model as the total mass profile, respectively. The 
vertical dotted line is the same as in Fig. 4. 
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Fig. 10. — The velocity dispersion profiles shown in Fig. 9 are converted to "temperature" 
via /cTdm = /^'^pC'DM- The ICM temperature profile derived in Sect. 4.3 and shown in 
Fig. 2 is overlaid with diamonds. "Temperatures" of member galaxies derived by firripag^i 
are also indicated with crosses, where the galaxy velocity dispersions (cxgai) are taken from 
observations by den Hartog & Katgert (1996). The star symbol with the gray vertical error 
bar indicates the central stellar velocity dispersion of the cD galaxy being converted to 
"temperature" by /^^pO'steiiar) where we use cisteiiar = 297 ± 12 (km s^^) as measured by 
Blakeslee & Tonry (1992). 



-al- 



io 



Q 

CO. 



10 




_l I I I I I I L. 



100 



Radius[arcsec] 



1000 




_] I I L I I I L. 



1000 



Fig. 11. — Radial profile of /?dm- The gray region and the hatched region with oblique lines 
correspond to the solutions shown in Fig. 9. The vertical dotted line is the same as in Fig. 4. 
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Fig. 12. — The hatched region with obhque hnes shows the velocity dispersion profile of the 
dark matter, when the NFW model is used as the total mass profile and the anisotropy of 
the velocity distribution is introduced as A = O.Gbjj^^^y^- The solution in the isotropic 
case is overlaid with the gray hatched region. The vertical dotted line is the same as in 
Fig. 4. 



- 33 - 




Radius[arcsec] 

Fig. 13. — Radial profile of AE. The gray region and the hatched region with oblique lines 
correspond to the solutions shown in Fig. 9. The vertical dotted line is the same as in Fig. 4. 
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Fig. 14. — (Left panel) Radiative cooling time of the ICM calculated from the density and 
temperature at each radius. (Right panel) Bolometric luminosity integrated within each 
radius. 
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Table 1. Best-fit parameters of tlie double (3 model fitted to the 0.8-10 keV brightness 

profile 



parameter 


value"^ 




^0,1 


(1.42 ±0.04) X 10-6 


[c/s/arcsec^] 




26.3 ± 1.0 


[arcsec] 


Pi 


0.562 ± 0.015 




5^0,2 


(1.92 ±0.12) X 10-^ 


[c/s/arcsec^] 


Rc,2 


106 ± 4 


[arcsec] 




0.633 ± 0.010 






83.2/74 





^Errors are 90% (A^^ = 2.7) confidence. 



Table 2. 



mass 
model 


P(),i 

[ Mq arcsec ^'^ ] 


Rc,i 

[arcsec] 


P(),2 

[ M0arcsec~"^ ] 


Rc,2 

[arcsec] 


ng{0) 
[cm~^] 




[arcsec (Mpc)] 


[10"Mq] 


1 King 

2 King 


1.28 X 10^ 
S.Slto lt X 10'' 


94.6 


2.74tJ;f| X 10« 


194.4111^, 


0.063 
0.086troS^ 


1124/77 
88.3/75 


1716t?i (2.041°;°^) 


t- OO+0.62 

J.OO_g gg 


NFW 


PO 

1.64t°;J| X 10^ 








0.1181°;°?^ 


88.7/77 


1958111 {2.33lO;°|) 


a ni+0-32 
«-U-'-_o.54 



*Mvir s |7r-Rj;j.po(l + z)^A, where po is the mean density of the universe at present, and A is the ratio between a cluster mean density 
and the mean density of the universe at the cluster-formation redshift z (=0.0616 for A1795) (see e.g. Kitayama & Suto 1996). 



